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Abstract 

We present adaptive finite element simulations of dendritic and eutectic solidification in binary and ternary alloys. 
The computations are based on a recently formulated phase-field model that is especially appropriate for modelling 
non-isothermal solidification in multicomponent multiphase systems. In this approach, a set of governing equations 
for the phase-field variables, for the concentrations of the alloy components and for the temperature has to be 
solved numerically, ensuring local entropy production and the conservation of mass and inner energy. To efficiently 
perform numerical simulations, we developed a numerical scheme to solve the governing equations using a finite 
element method on an adaptive non-uniform mesh with highest resolution in the regions of the phase boundaries. 
Simulation results of the solidification in ternary Ni6oCu4o_ :E Cr 3; alloys are presented investigating the influence of 
the alloy composition on the growth morphology and on the growth velocity. A morphology diagram is obtained 
that shows a transition from a dendritic to a globular structure with increasing Cr concentrations. Furthermore, we 
comment on 2D and 3D simulations of binary eutectic phase transformations. Regular oscillatory growth structures 
are observed combined with a topological change of the matrix phase in 3D. An outlook for the application of our 
methods to describe AlCu eutectics is given. 
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1. Introduction transformations, complex thermodynamic interac- 

tions and pattern formations. The complex phase 
Multicomponent alloys form the most important diagrams of multicomponent alloys show a vari- 

class of metallic materials for technical and in- et y of different solidification processes such as eu- 

dustrial processes. Combined with the number of tectic > peritectic and monotectic types of reaction. 

components is a wealth of different phases, phase Since the most common solidification microstruc- 

tures occuring in industrial alloys, such as Al and 
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gations on these two types of structures. 

The phase-field modelling technique has made a 
significant progress during the last decade in simu- 
lating complex microstructures, such as dendritic 
and eutectic growth; see e.g. [1,2]. The newly intro- 
duced variable 4>(x, t) in these models indicates the 
physical state of a system at each point in space. 
Differing from sharp interface concepts, </>(x, t) has 
a smooth transition at the interfaces (phase bound- 
aries). In the bulk phases, <j>(x, t) takes on constant 
values, e.g. 4>(x, t) = for the liquid and 4>(x, t) = 
1 for the solid phase. The governing equations con- 
sist of parabolic partial differential equations for 
the evolution of the phase state coupled to mass 
and heat transport equations. They serve for per- 
forming simulations of complex growth morpholo- 
gies without explicitly tracking the phase bound- 
aries. In a recent article [3] , the phase-field method- 
ology has been generalized to the case of arbitrary 
numbers of N components and M phases in alloy 
systems by introducing a vector of concentrations 
c = (ci, . . . , cat) and a vector of phase-field vari- 
ables (j) = (0i, . . . , (j> M ). 

For computations based on diffuse interface 
models, it is required that the spatial resolution 
of the numerical method must be greater than the 
thickness of the diffusive phase boundary layer. 
The interfacial thickness itself must be less than 
the characteristic scale of the growing microstruc- 
turc. In this case, a non-uniform grid with adaptive 
refinement can dramatically reduce the use of com- 
putational resources against a uniform grid with 
the same spatial resolution. For these reasons, we 
use an adaptive finite element method for solving 
the governing equations of the new multicompo- 
nent phase-field model. 

The present report is part of our intension 
to apply the phase-field model and the numer- 
ical method to simulate solidification processes 
in real alloys which are usually multicomponcnt. 
We present simulation results of ternary dendritic 
and of binary eutectic growth. In particular, the 
ternary Ni6oCu4o- x Cr x alloy system is considered 
to investigate the influence of interplaying solute 
fields on the interface stability, on the growth veloc- 
ity and on the characteristic type of morphology. 
By choosing the Ni-Cu-Cr system, we build upon 
the binary Ni-Cu system which has been explored 



by phase-field modelling (e.g. [4]) and by molec- 
ular dynamics simulations in several papers (e.g. 
[5] ) . Hence, physical parameters are relatively well 
established. For our application to eutectic solidi- 
fication, we constructed a symmetric binary model 
phase diagram. Special emphasis is laid on a dis- 
cussion of oscillatory growth structures. 



2. Phase-field model 

The phase-field model in [3] has been formulated 
in a thcrmodynamically consistent way and allows 
for an arbitrary number of phases and components. 
It is defined solely via bulk free energies of the in- 
dividual phases, the surface energy densities of the 
interfaces, the diffusion and mobility coefficients 
and yields classical moving boundary problems in 
the sharp interface limit. The approach is based on 
an entropy density functional S(e, c, <fi) of the form 



s(e, c, (j)) - ( ea{4>, V0) + -w((f>) 



dx. 
(1) 

The bulk entropy density s depends on the inter- 
nal energy e, on the concentrations of components 
c = (ci)Hi, and on the phase-field variable <p = 
(</)a)a=i- The thermodynamics of the interfaces is 
given by the second and third term in the integral 
of Eq. (1) and is determined by the gradient energy 
density a(</>, V0), the multi well potential w((f>) and 
a small scale parameter e related to the thickness 
of the interface. The gradient energy and the multi 
well potential depend on the surface energy density 
a and its anisotropy can be taken into account us- 
ing appropriate choices of a(<j), V</>). The variable 
<p a with < <p a < 1 denotes the local fraction of 
phase a. It is required that the concentrations and 
phase-field variables fulfill the constraints 



N 



M 



(2) 



i=i 



The evolution equations for the phase fields are 
postulated as 

1 / 

u>ed t (t> a = e(Va V Q - a </>J - -w,<t> a - LJ ^r - \ 

(3) 
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where f(T, c, <f)) is the free energy density and A 
is a Lagrange multiplier such that the constraint 
in Eq. (2) for phase fields is satisfied. The kinetic 
factors = u)((j), V0) describes anisotropic interface 
kinetics and is related to the kinetic coefficient [i of 
atomic attachment in the linear response function 
"growth velocity-interfacial undercooling" at a flat 
front. 

Considering an ideal solution system, we define 
the free energy density as follows 

N M _ „„ N 

f(c, d>) = ]T ]T Ci L? — ^h(cj> a )+— £ Cl In a, 

Ha=l ±i Vm i=l 

(4) 

where Lf and T" are the latent heat and the melt- 
ing temperature of the component i in the phase a, 
respectively. R is the gas constant, v m is the mo- 
lar volume and h((f> a ) is a monotone function on 
the interval [0, 1] satisfying h(0) = and h(l) = 1. 
Additionally, we consider an isothermal approxi- 
mation with the temperature T being constant. 

Assuming that the mass fluxes are linear func- 
tions of the thermodynamic driving forces, i.e. 
chemical potentials in isothermal approximation, 
mass balance equations can be written as 

d t c i = -v(^L ij (c,<t>)V=^j , (5) 

with chemical potentials in — / ;Cj and mobility 
coefficients given by 

U (c, 4>) = yAc, ( <% - — ] • (6) 

The form of Eq. (6) allows different values of the 
bare trace diffusion coefficients Di{(f>) for the dif- 
ferent components i and satisfies the constraint for 
concentrations in Eq. (2). 

The evolution Eqns. (3) and (5) are solved us- 
ing a finite element method with Lagrange ele- 
ments and linear test functions. For the time evo- 
lution, a semi-implicit formulation is discretized 
to achieve better numerical stability. We generate 
a non-uniform adaptive mesh having the highest 
order of spatial resolution in the vicinity of the 
solid-liquid interface where the gradients of phase 
fields and concentrations reach maximal values. 
The mesh structure is adopted in time according to 



the evolution of phase and concentration fields us- 
ing a refinement criterion based on the phase-field 
and concentration gradients. 



3. Dendritic growth in ternary alloys 

By numerical simulations, we investigate how 
the change of the alloy composition influences the 
growth morphology for a given solidification con- 
dition. In particular, we fix the initial undercool- 
ing and choose a ternary Ni 6 oCu4o- £! ;Cr a ; alloy as 
a prototype for this study. To recover the corre- 
sponding solidus and liquidus lines of the binary 
Ni-Cu and Ni-Cr phase diagrams in the region of 
concentrations up to 40at.% of Cu or Cr, the fol- 
lowing physical parameters are used: Melting_tem- 
peratures T m = 1728 K, T Cu = f358 K, T Cr = 
1465 K; latent heats L m = 2350 J/cm 3 , L Cu = 
1728 J/cm 3 , Z Cr = 1493 J/cm 3 and a molar vol- 
ume v m = 7.42 cm 3 . Since we do not consider the 
complete eutectic Ni-Cr phase diagram, the values 
of the melting temperature and of the latent heat 
for Cr are not the real physical data, but adjustable 
parameters in order to recover the actual binary 
phase diagram in the given region of concentra- 
tions. They are marked by a tilde to emphasize this 
difference from the data for Ni and Cu. Applying 
these values leads to a partition coefficient k e = 
0.843, to a liquidus slope m e = —3.27 K/at.%, and 
to a freezing range AT = m e (k e ~ l)/fc e Coo = 
24 K for binary Ni 60 Cu4o. Similarly, we obtain k e — 
0.905, m e = -2.08 K/at.%, and AT = 8.7 K for 
the binary Ni 60 Cr4 system and the corresponding 
equilibrium phase diagram. 

We assume that both surface properties, the sur- 
face energy density a and the kinetic coefficient \i 
do not depend on the alloy composition and have 
the values a = 0.37 J/m 2 and [i = 3.3 mm/(s 
K), refcring to [4]. The anisotropy of the interface 
properties plays an important role in the selection 
of the operating state during dendritic growth. In 
this study, we use the values calculated from molec- 
ular dynamics simulations. In [5], the strength of 
the surface energy anisotropy is given as 0.023 and 
the strength for kinetic anisotropy as 0.169. 

The diffusion coefficients in the melt are Z?Ni = 
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3.82 x 1CT 9 m 2 /s [5], D Cu = 3.32 x 10~ 9 m 2 /s 
[5] and Z?o = 1-5 x 10~ 9 m 2 /s. The diffusion co- 
efficients in the solid phase are set as an equal 
value 1CP 13 m 2 /s for all components. The small 
length scale parameter e in the entropy functional 
(Eq. (1)) is e = 0.1 /urn. 

Using the physical parameters given above, 
we carried out a series of numerical computa- 
tions for different alloy compositions varying from 
Ni6oCu36Cr4 to Ni6oCu4Cr36. The concentration 
of Ni was kept at 60 at.% and the initial undercool- 
ing was fixed at 20 K measured from the equilib- 
rium liquidus line in the phase diagram at a given 
composition of the melt. A morphological transi- 
tion from dendritic to globular growth occurs at a 
melt composition of about NigoCu2oCr2o- The left 
side of Fig. 1 shows the dendritic morphologies ob- 
served for Cr concentrations less than 20 at.%. The 
right side of the figure displays globular morpholo- 
gies for Cr concentration crossing this threshold. 
The velocity of the dendritic/globular tip increases 
linearly from 1.19 cm/s to 3.24 cm/s with increas- 
ing the concentration of Cr. An analogous type 
of dendritic to globular morphology transition de- 
pending on the undercooling has been observed in 
numerical simulations of binary alloy solidification 
and is discussed in [6]. 

Qualitatively such transitions can be expected 
by comparing the initial undercooling AT with the 
freezing range AT : For the binary Ni 6 oCu4o al- 
loy the undercooling AT = 20 K is less than the 
freezing range ATo = 24 K, whereas for the binary 
Ni 6 oCr 4 o system the undercooling AT = 20 K is 
greater than ATo = 8.7 K. Thus the morphology 
change is related to the transition from a two-phase 
region (above the solidus line) to a one-phase re- 
gion (below the solidus) in correspondance with 
the phase diagram. Qualitative analysis of growth 
structures should take into account the kinetic ef- 
fects of the solid-liquid interface as well as the sur- 
face tension acting as stabilizing factor and the dif- 
fusion properties of alloy components as destabi- 
lizing forces due to concentration gradients [7] . 



4. Simulation of binary eutectic growth in 
2D and 3D 

To simulate eutectic phase transitions in a bi- 
nary A-B alloy, where two solid phases a and 
j3 grow into an undercooled melt, we have con- 
structed a typical eutectic phase diagram via the 
method of common tangents. In our computations, 
we consider the widely observed phenomenon of 
regular oscillations along the solid-solid interface 
driven by the triple junctions (Fig. 2). The initial 
concentration of the melt is chosen at the eutectic 
composition ca = cb = 0.5, where the two solid 
phases grow with equal phase fractions. A transi- 
tion from steady state lamellar growth to an oscil- 
latory pattern formation is found for an increas- 
ing difference of the initially set up phase fraction. 
For a ratio of 1:3 between the initial phase frac- 
tion of solid a:f3, lamellar growth is re-established. 
For ratios of 1:4 and 1:6, regular oscillations with 
a characteristic amplitude and wave length as in 
Fig. (2) can be observed, whereas for a ratio 1:7, 
the /? phase overgrows the initially dominating a 



Ni 6 oCu 28 Cr 12 



Ni 60 Cu 12 Cr 28 




Ni 60 Cu 24 Cr 16 



Ni 60 Cu 16 Cr 24 



Fig. 1. Dendritic to globular morphological transition for 
different alloy compositions. The atomic percents of Cu 
and Cr are exchanged keeping Ni fixed at 60 at.%. The 
shaded regions correspond to the solid phase and the solid 
lines represent the isolincs of average concentration of Ni 
in the solid phase. 
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Fig. 2. Regular oscillations along the solid-solid interface 
during binary eutectic growth driven by the motion of the 
triple junction in 2D. 




Fig. 3. Topological change of the microstructure with a 
rods embedded in a /3 matrix phase and visa versa. The 
formation results from regular 3D oscillations along the 
solid-solid interface. 

solid phase. In such nucleation event takes 

place. Two-dimensional oscillatory structures have 
been discussed in [8] by boundary integral method 
and more recently in [9] by means of phase-field 
modelling. 

The analogous type of oscillation can be ob- 
served for eutectic microstructure formations in 
three dimensions, Fig. 3. Performing an alternat- 
ing topological change, a solid rods are embed- 
ded in a (3 matrix followed by the opposite situa- 
tion of /3 crystals embedded in an a matrix. Fur- 
ther three-dimensional simulations of eutectic mi- 
crostructures are reported in [10,11] 



5. Conclusion 

A new phase-field model for non-isothermal so- 
lidification in multicomponent multiphase alloy 



systems has been applied to ternary Ni-Cu-Cr pri- 
mary phase growth and to binary eutectic growth 
of two solid phases into an undercooled melt. A 
transition from a dendritic to a globular structure 
is observed in our simulation series, while succes- 
sively exchanging Cu by Cr. The experiences with 
a binary model eutectic will enable us to apply the 
modelling and simulation technique in eutectic and 
off-eutectic Al-Cu alloys in a forthcoming paper. 
Steady-state lamellar domains as well as different 
types of oscillations persisting over a wide range of 
compositions and growth rates are reported in [12]. 
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